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Abstract 

There are a number of constraints which limit the current and voltages which 
can be applied on a multiple drive electrical imaging system. One obvious con- 
straint is to limit the maximum Ohmic power dissipated in the body. Current 
patterns optimising distinguishability with respect to this constraint are singular 
functions of the difference of transconductance matrices with respect to the power 
norm, (the optimal currents of Isaacson). If one constrains the total current (L 1 
norm) the optimal patterns are pair drives. On the other hand if one constrains the 
maximum current on each drive electrode (an L°° norm), the optimal patterns have 
each drive channel set to the maximum source or sink current value. In this pa- 
per we consider appropriate safety constraints and discuss how to find the optimal 
current patterns with those constraints. 



1 Introduction 

The problem of optimizing the drive patterns in EIT was first considered by Seagar [jl|] 
who calculated the optimal placing of a pair of point drive electrodes on a disk to max- 
imize the voltage differences between the measurement of a homogeneous background 
and an offset circular anomaly. Isaacson [[|], and Gisser, Isaacson and Newell [Q] 
argued that one should maximize the L 2 norm of the voltage difference between the 
measured and calculated voltages constraining the L 2 norm of the current patterns in 
a multiple drive system. Later [Q] they used a constraint on the maximum dissipated 
power in the test object. Eyoboglu and Pilkington [0| argued that medical safety legis- 
lation demanded that one restrict the maximum total current entering the body, and if 
this constraint was used the distinguishability is maximized by pair drives. Cheney and 
Isaacson [||| study a concentric anomaly in a disk, using the 'gap' model for electrodes. 
They compare trigonometric, Walsh, and opposite and adjacent pair drives for this case 
giving the dissipated power as well as the L 2 and power distinguishabilies. Koksal and 
Eyoboglu [||] investigate the concentric and offset anomaly in a disk using continuum 
currents. 

Yet another approach [^J is to find a current pattern maximizing the voltage differ- 
ence for a single differential voltage measurement. 

* Department of Mathematics. UMIST, UK 

^ Department of Applied Physics, University of Kuopio, Finland 

* School of Engineering, Oxford Brookes University.UK 



1 



2 Medical Electrical Safety Regulations 



We will review the current safety regulations here, but notice that they were not de- 
signed with multiple drive EIT systems in mind and we hope to stimulate a debate 
about what would be appropriate safety standards. 

For the purposes of this discussion the equipment current ("Earth Leakage Current" 
and "Enclosure Leakage Current") will be ignored as the emphasis is on the patient 
currents. These will be assessed with the assumption that the equipment has been 
designed such that the applied parts, that is the electronic circuits and connections 
which are attached to the patient for the delivery of current and the measurement of 
voltage, are fully isolated from the protective earth (at least 50MQ.). 

IEC601 and the equivalent BS5724 specify a safe limit of 100 /jA for current flow 
to protective earth ("Patient Leakage Current") through electrodes attached to the skin 
surface (Type BF) of patients under normal conditions. This is designed to ensure that 
the equipment will not put the patient at risk even when malfunctioning. The standards 
also specify that the equipment should allow a return path to protective earth for less 
than 5 mA if some other equipment attached to the patient malfunctions and applies 
full mains voltage to the patient. Lower limits of 10 /jA (normal) and 50 /jA (mains 
applied to the patient) are set for internal connections, particularly to the heart (Type 
CF), but that is not at present an issue for EIT researchers. 

The currents used in EIT flow between electrodes and are described in the standards 
as "Patient Auxiliary Currents" (PAC). The limit for any PAC is a function of frequency, 
100 microamps from 0.1 Hz to 1 kHz; then 100/ /jA from 1 kHz to 100 kHz where / 
is the frequency in kHz; then 10 mA above 100 kHz. The testing conditions for PAC 
cover 4 configurations; the worst case of each should be examined. 

1. Normal conditions. The design of single or multiple current source tomographs 
should ensure that each current source is unable to apply more than the maximum 
values given. 

2. The PAC should be measured between any single connection and all the other 
connections tied together, a) if the tomograph uses a single current source then the 
situation is similar to normal conditions (above) b) if the tomograph uses multiple 
current sources then as far as the patient is concerned the situation is the same as normal 
conditions. The design of the sources should be such that they will not be harmed by 
this test. 

3. The PAC should be measured when one or more electrodes are disconnected 
from the patient. This raises issues for multiple-source tomographs : a) if an isolated- 
earth electrode is used then the current in it will be the sum of the currents which should 
have flowed in the disconnected electrodes; they could all be of the same polarity. The 
isolated-earth electrode should therefore include an over-current sensing circuit which 
will turn down/off all the current sources, b) If no isolated-earth electrode is used then 
the situation is similar to normal conditions. 

4. The PAC should be measured when the disconnected electrodes are connected 
to protective earth. This introduces no new constraints given the tomograph is fully 
isolated. 

3 Constrained Optimization 

Let V = (Vi , . . . , Vk ) t be the vector of potentials measured on electrodes when a pat- 
tern of currents / = (I\ , . . . ,Ik ) T is applied. These are related linearly by R the transfer 
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impedance matrix: V = RI. For simplicity we will assume the same system of elec- 
trodes is used for current injection and voltage measurement. We will also assume that 
the conductivity is real and the currents in-phase to simplify the exposition. A model 
of the body is used with our present best estimate for the conductivity and from this we 
calculate voltages V c for the same current pattern. Our aim is to maximize the distin- 
guishability ||V— V c 1 1 2 = ||(^ — -Rc)/[[2- The use oftheL 2 norm here corresponds to the 
assumption that the noise on each measurement channel is independent and identically 
distributed. If there were no constraints on the currents the distinguishability would be 
unbounded. 

The simplest idea [||] is to maximize \\(R — R^)I\\2 subject to ||/||2 < M for some 
fixed value of M. The solution of this problem that / is the eigenvector of R — R c corre- 
sponding to the largest (in absolute value) eigenvalue. One problem is that the 2-norm 
of the current has no particular physical meaning. In a later paper [|J it was proposed 
that the dissipated power be constrained, that is / • V = I 1 RI. The optimal current is 
the eigenvector of (R — R c )R~ l l 2 . (The inverse implied in the expression R~ x l 2 has to 
be understood in the generalized sense, that is one projects on to the space orthogonal 
to (1,...,1) and then calculates the matrix exponent 1 /2.) In practical situations in 
medical EIT the total dissipated power is unlikely to be an active constraint, although 
local heating effects in areas of high current density may be an issue. Even in indus- 
trial applications of EIT, the limitations of voltages and currents handled by normal 
electronic devices mean that one is unlikely to see total power as a constraint. One 
exception might be in EIT applied to very small objects. 

As we have seen a reasonable interpretation of the safety regulations is to limit 
the current on each electrode to some safe level 7 max . We will refer to this as an L°° 
constraint. This corresponds to a convex system of linear constraints — 7 max < 4 < /max- 
When we maximize the square of the distinguishabilty, which is a positive definite 
quadratic function of /, with respect to this set of constraints it can be seen that the 
maximum must be a vertex of the convex poly tope {/ : maxfc{|4|} = /max, Lit 4 = 0}. 
For example, for an even number 2« of electrodes the 2 "C„ vertices are the currents with 
each 4 = ±/ lnax , and an equal number with each sign. For the circularly symmetric case 
these are the Walsh patterns referred to in [|]| . 

If one wanted to be safe under the multiple fault condition that all the electrodes 
driving a current with the same sign became disconnected, and the safety mechanism 
on the isolated-earth failed, one would employ the L 1 constraint |4| < 2/ ma x- Again 
this gives a convex feasible set. In this case a polyhedron with vertices / such that all 
but two 4 are zero, and those two are / ma x and — / ma x- These are the pair drives as 
considered by Seagar. Pair drives were also considered by [§|,[0],[0] for single circular 
anomalies. Notice that L 1 optimal currents will be pair drives for any two- or three- 
dimensional geometry and any conductivity distribution. 

Another constraint which may be important in practice is that the current sources 
are only able to deliver a certain maximum voltage V m ax close to their power supply 
voltage. If the EIT system is connected to a body with transfer impedance within its 
design specification then the constraints —Vmax < Vk < Vmax will not be active. If they 
do become active then the additional linear constraints in / space — V m ax <R~ l I < V max 
(here R~ l is to be interpreted as the generalized inverse), will still result in a convex 
feasible region. 

When any of the linear constraints are combined with quadratic constraints such as 
maximum power dissipation the feasible set of currents is still convex but its surface is 
no longer a poly tope. 



3 





Single anomaly 


Two anomalies 


Constraint 


Voltage diff. 


Power 


Voltage diff. 


Power 


L l Best pair drive 
1? optimal 
Power optimal 
L°° optimal 


440.0911 
812.3243 
518.9126 
1653.2673 


1200.7618 
3356.3035 
1768.3048 
19261.9208 


347.3579 
571.9161 
321.616 
968.2656 


1185.4935 
3170.0985 
1352.6896 
16798.3363 



Table 1 : L 2 norm of voltage difference, and dissipated power for one and two anoma- 
lies with a variety of constraints. The constraint levels have been chosen so that the 
maximum electrode current is the same on each 

4 Numerical Results 

Although we can easily find the vertices of the feasible region there are too many for 
it to be wise to search exhaustively for a maximum of the distinguishability. For 32 
electrodes for example there are 32 Ci6 > 6 x 10 8 . Instead we use a discrete steepest 
ascent search method of the feasible vertices. That is from a given vertex we calculate 
the objective function for all vertices obtained by changing a pair of signs, and move to 
whichever vertex has the greatest value of the objective function. For comparison we 
also calculated the L 2 optimal currents, the optimal currents for the power constraint, 
and the optimal pair drive (L 1 optimal). 

We used a circular disk for the forward problem, and the EIDORS Matlab tool- 
box [^] for mesh generation and forward solution. The mesh and conductivity targets 
can be seen in Figure @. Our results are interesting in that for the cases we have stud- 
ied so far the L°° optimal currents have only two sign changes. The distinguishabilies 
given in Table [l] should be read with caution, as it is somewhat unfair to compare for 
example power constrained with L°° patterns. They are designed to optimise different 
criteria. However the contrast between pair drive and L°° is worth noting as the majority 
of existing EIT systems can only drive pairs of electrodes. 

The greatest current densities occur at the contact points between the electrode 
boundaries and skin. At each electrode this current density is determined mainly by 
the total electrode current, the contact impedance and the skin conductivity just below 
the electrode. These factors dominate the current density near the electrode boundaries 
and the other electrode's currents have a much smaller contribution to the maximum 
current densities 

5 Conclusions 

If using optimal current patterns one should be sure to use the right constraints. We 
suggest that in many situations the L°° constraint may be the correct one. We have 
demonstrated that it is simple to compute these optimal patterns, and the instrumenta- 
tion required to apply these patterns is much simpler than the L 2 or power norm pat- 
terns. While still requiring multiple current sources, they need only be able to switch 
between sinking and sourcing the same current. 
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Figure 3: Mesh and conductivity anomalies. 
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